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ABSTRACT 

Primordial magnetic fields generated in the early universe are subject of consider- 
able investigation, and observational limits on their strength are required to constrain 
the theory. Due to their impact on the reionization process, the strength of primor- 
dial hclds can be limited using the latest data on reionization and the observed UV- 
luminosity function of high-redshift galaxies. Given the steep faint-end slope of the 
luminosity function, faint galaxies contribute substantial ionizing photons, and the 
low-luminosity cutoff has an impact on the total budget thereof. Magnetic pressure 
from primordial fields affects such cutoff by preventing collapse in halos with mass 
below 10 10 M Q (i?o/3 nG) 3 , with Bo the co- moving held strength. In this letter, the 
implications of these effects are consistently incorporated in a simplified model for 
reionization, and the uncertainties due to the cosmological parameters, the reioniza- 
tion parameters and the observed UV luminosity function are addressed. We show 
that the observed ionization degree at z ~ 7 leads to the strongest upper limit of 
Bq<>2 — 3 nG. Stronger limits could follow from measurements of high ionization de- 
gree at z > 7. 

Key words: magnetic helds - galaxies: high-redshift - intergalactic medium - cos- 
mology: theory - dark ages, reionization, first stars 



1 INTRODUCTION 



The origin of magnetic fields in the universe is subject of ac- 
' tive study. Magnetic f ields are a comm o n feature of nearb 
& \ and distant g alaxies (|Beck et al.1 Il99fj ; iBernet et ail |200 
iMurphvl 12009). where the y quickly build u p due to small 



201ll: ISchlickeiser fc Shuklal |2003l: iBertone et al l 120061 : 
Ando et al.ll20ld ). In particular JMiniati fe Belli (|201lh make 



scale dynamo action (e.g 



l20ld : iFederrath et al J 12011 



Schleicher et al 



Schober et al 



l20ld : ISur et al.1 
20111 ). Magnetic 



fields exist also in the intergalactic medi um (IGM). In par - 
ticular, galaxy clusters exhibit fj,G fields ( Clarke et al.ll200ll) 



consistent predictions for the magnetic fields observed 
in the co smic voids. Primordial mo dels of magnetogen- 
esis (e.g. iGrasso fc Rubinstein! l200ll l provide an alterna- 
tive scenario. Unlike most astrophysical cases, their seeds 
are not-necessarily negligibl e. This is particularly so for 
some inflationary scenarios dTurner fc Widrowl [l988) , elec- 
troweak ph ase transition jBavm et al.lll996 | ). or QCD phase 



which require non-negl i gible initial seeds | Dubois fc Teyssien 
120081 : iRvu et al.1 120081 : iMiniati fc Martini l201ll )~ 



troweak pnasc transition (rjaym ct al. 1990), or you pliase 
transition dQuashnock et al.l Il989l; I Cheng fc Olintol Il994l : 



but 



Also, tiny 

magnetic fields appear to exist in cos- 
as recently suggeste d by gamma-ra y ex- 



significant 
voids 

periments iTNeronov fc Vovkll2oTol: iTavecchio et al.l l20ld : 

The gen- 



iTavlor et alj|2011 



but see 



Broderick et al. 



(|2Q11|)). 



eration of magnetic fields in cosmic environments is gener- 
ally slower than in galaxies, owing to the longer dynamical 
timescales. 

Various astrophysical mechanisms have been proposed 
for the origin of intergalactic magnetic fields (jMiniati fc Belli 
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ISiel et al.l[l997l : iBaneriee fc Jedamzildl2003h . Therefore 
eral efforts to constrain primordial seeds have been made. 

Probing magnetic fields at high redshift is however 
still challenging. A primordial magne tic field has an im- 
pact on the CMB anisotropies (e.g. iBarrow et al.l 1 19971 : 
ISubramanianl l200d : iDurrerl l2007u and the latest analysis 
implies an upper limi t on the co- moving field strength 
B <3nG at 95% CL l|Yamazaki et alj|2010h . where B = 
B(z)/(1 + z) 2 . Constraints from Big Bang nucleosynthesis 
can a lso be obtained, with B^lfiG l|Grasso fc Rubinstein! 
1996). Considerations of gravitational wave production 
seem to rule out magnet ic field production during pri- 
mordial phase transitions (|Caprini et a l. 2009). Finally, as 
pointed out in various works (e.g. IColes! Il992l : | Kim et al.l 
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19961; Battancr ot al. I Il997l I Seth i fc SubramaniarJ |2005; 
Tashiro fc Sugivama 20061 ), the magnetic field can affect the 



evolution of the IGM and the growth of structure. These 
effects can be used to set indepen dent constraints on pr i- 
mordial fields during reionization l|Schleicher et al.l t2008). 
In particular, the magnetic field adds a pressure term which 
affects the Jeans mass. Also, if the magnetic field is not par- 
allel to the electric current, charges are subject to a Lorentz 
force. This gives rise to ambipolar diffusion of the charges 
through the neutrals and heating of the IGM due to colli- 
sions between the two species. 

In this letter, we constrain primordial magnetic fields 
by considering their impact on the reionizaton process. Our 
analysis is based on the latest data of the high redshift IGM. 
In particular, we use observed UV lumino s ity fu nction at 
redshifts z ~ 4 - 8 dBouwens et al.1 12010| . 1 20 111) and the 
recen t observations of a z — 7.085 quasar (Mortlock et al. 
l201ll ) which indicates a fraction of neutr al hydrogen of 
10" 4 - 10" 3 at z ~ 7 (|Bolton et al.l l201lT l. We find that 
Bp <^2 — 3nG, similar to fi ndings from the latest CMB analy- 
sis l|Yamazaki et alj|2010l ). More importantly, our constraint 
is complementary with the latter as it is based on differ- 
ent physical processes and refers to a different cosmological 
epoch. Th roughout we use the c osmological parameters from 
WMAP7 (|Komatsu et al.ll201ll ). 



2 TECHNICAL APPROACH 

In order to assess the impact of primordial magnetic 
fields on reionization we solve the following equation for 
the ev olution of the volu me fraction of ionized hydrogen, 
Qhii l|Madau et al.lll999D 



dQun 
dt 



t r . 



mi SFR(z)/, 



nff(0) 



(1) 



where uh{0) is the comoving number density of neutral hy- 
drogen and the other parameters are described below. For 
ou r reference model we assume the parameter values as 
in lBouwens et ail (|201ll . see also references therein), which 
reproduce the Thomson optical depth observed by WMAP7 
within la. However, for the cosmo logical parameters w e 
use the more recent WMAP7 results (|Komatsu et al.ll201ll ). 
As discussed below, this makes a negligible difference. So, 
the adopted escape fraction of Lyman-continuum photons 
is /esc = 0.2. The temperature of the reionized IGM is 
2 x 10 4 K, leading to an hydrogen recombination timescale 



1.0 Gyr 



7 



(2) 



where C3 is the clumping factor normalized to a value of 
3. Note that the recombination timescale would decrease 
by a factor of 2 for an IGM temperature of 10 4 K. The 
factor 10 53 ' 2 in Eq. (fjTJ is the production rate of Lyman- 
continuum photons per Moyr -1 of forming stars, enhanced 
by 30 % for metal-poor st ars in galaxies in the early uni- 
verse (jBouwens et al.f201ll ). The cosmic star formation rate, 
SFR, uncorrected for dust attenuation and in Moyr~\ is 



SFR 



)(Muv)f>*(Muv)d (iQ-^Ww-M^ ^ (3) 



where &(Muv) is the Schechter function as a function of 
the UV magnitude, Muv, and M V v its characteristic mag- 
nitude. As for the parameters entering the Schechter func- 
tion (inclu ding M* rv ) we use t he observational values in 
Table 1 of[g ouwens et al.1 l|201ll ). The star formation rate 
(uncorrected for dust attenuation) for a galaxy with mag- 
ni tude Muy is o btain ed by manipulation of the expression 
in lMadau et al.l (1 19981 ) and reads 



5.43 x 10 



-(M uv +20)/2.5 



M yr \ 



(4) 



Eq. Q contains contr i bution s up to a limiting magnitude, 
M^y. lBouwens et al.1 <|201ll ) assumes M^ x = -10, corre- 
sponding to the mass scales where galaxy formation is sup- 
pressed due to inefficient gas cooling and/or feedback effects. 
Here we generalize this limit to include effects of magnetic 
pressure on the Jeans mass. For this purpose we use the fil - 
tering mass M F (|Gnedin fc Huil 19981 ; [Schleicher et al.ll2008l ) 



M 



2/3 



da'M 2 g /3 {a) 



1/2 



(•5) 



where, M g denotes the maximum of the thermal, Mj, and 
magnetic Jeans masses, Mf. These are, respectively 



Mj = 2Ma 



0.2 km/s i V 1000 cm 



(6) 



with eg the sound speed and n t he gas number density, 
and (|Subramanian fc Barrowl ll998) 



Mf = 10 



1O M ( 



Bp 
3 nG 



(V 



with Bo the co-moving field strength. Finally, to convert the 
filtering mass scale into a UV magni tude, we use the f ollow - 
ing relation inferred from Fig. 10 of ISalvaterra et al] (|201ll ) 
(which also consistently reproduces the above limiting mag- 
nitude in absence of magnetic field) 

M uv (M h&lo ) = -20 - 3 x log 10 . (8) 

In addition to the volume fraction of ionized hydrogen, 
our calculation follows the time evolution of the temperature 
and ionization degree of the partially-ionized IGM, and the 
magnetic field strength . The full set of equations is given in 
ISchleicher et al.l J2008I) an d was implemented in the REC- 
FAST coded | Seager et al.lll999l ). In particular, the magnetic 
energy decreases with time due to adiabatic expansion and 
dissipation into heat through ambipolar diffusion. For strong 
fields this can change substantially the thermal and ioniza- 
tion history of the IGM, although full reionization cannot 
be achieved this way, as collisions induced by ambipolar dif- 
fusion become inefficient at high temperatures. 

With the calculation results, the Thompson scattering 
optical depth is computed as 



n H (0)c 
Ho 



(1 + 



x eff (z)a T - . 

y'^A + fim(l + zY 



dz, (9) 



where <jt is the Thompson scattering cross section and the 
effective ionization degree is 



Qhii + (1 - Qmi)xe 



(10) 



http: / /www. astro. ubc.ca/people/scott/recfast. html 
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Figure 1. The reionization optical depth as a function of the co- 
moving field strength. Upper panel: Variation of the cosmological 
parameters. Mid panel: Variation of the reionization parameters. 
Lower panel: Variations in the Schechter function. The curves 
used to derive the upper limit are marked with a thick line. 



Figure 2. The ionization degree at z = 7 as a function of the co- 
moving field strength. Upper panel: Variation of the cosmological 
parameters. Mid panel: Variation of the reionization parameters. 
Lower panel: Variations in the Schechter function. The curves 
used to derive the upper limit are marked with a thick line. 



contributed by both the fully ionized volume fraction as well 
as the partially ionized gas. 



3 RESULTS 

In the following, we set the upper limits on the comov- 
ing magnetic field strength, Bo, using the Thomson optical 
depth (a redshift-integrated quantity) as well as the the ob- 
served ionization degree at z ~ 7. We discuss the uncertain- 
ties due to the cosmological parameters, the reionization pa- 
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Figure 3. The ionization degree at different redshifts, as a func- 
tion of the co-moving field strength (in our reference scenario). 

rameters and the observed UV luminosity function. Finally, 
we discuss the possibility of more stringent constraints from 
future progress on the high redshift IGM. 

3.1 Constraints from the Thompson optical depth 

Fig. [1] shows the constraints on Bo due to the observed 
optical depth, r e , and its uncertainties. Starting from our 
reference case in the top panel, we see that for weak mag- 
netic fields the optical depth is well within the 1<t range of 
WMAP7 results. As Bo increases the low-luminosity cut- 
off of halos contributing ionizing photons is raised and at 
Bo^ nG, r e starts to drop. It reaches a minimum at several 
nG and then, as collisional ionization resulting from heating 
due to ambipolar diffusion kicks in, r e increases again (how- 
ever, values of several n G for Bp, would be inco nsistent with 
CMB observations, e.g., lYamazaki et al.ll2010h . This is the 
qualitative trend in all panels of Fig. 1, 2 and 3. 

The top panel also shows the effects of 3<r variations of 
the cosmological parameters, on r e . For the purpose of the 
current discussion, higher values of Ho and smaller values of 
Qm increase the value of r e and thus weaken the constrain 
on Bo. Thus, the constraint on Bo is set by the highest 
curve, which says that the required field strength at the 
95% CL limit is 3 nG. Not e also that adopting a cosmology 
as in lBouwens et all (|201ll ) would have minor effects on our 
results. 

In the mid panel we consider instead the effect due to 
uncertainties in the reionization parameters. Such uncertain- 
ties can in principle be very large. However, a large clumping 
factor ~ 30 is inconsistent with the observed ionization frac- 
tion (see Fig. 2) even in the limit of weak magnetic fields. 
Likewise, a very large escape fractio n / e3C ~ 50% produce 
too h i gh an ionization rat e at z<^Q (|Srbinovskv fc Wvithel 
|2010| ; iMiniati et all 120041 ). With these considerations, a 
meaningful 95% CL upper limit of 2nG is implied by our 
results. 

Finally, in the bottom panel we consider the uncertain- 
ties due to 3<J observational errors in the parameters a, M^ v 
and $*, of the Schechter function. Here, it is large values of 
$* or MJy that increase r e , and therefore weaken our con- 



straint. Thus, the highest curve on the plot shows that the 
computed value of r e remains within 2a of the WMAP7 
result for BoSS3nG. Note that the redshift evolution of the 
Schechter function parameters is also very uncertain, but we 
have verified that it does not lead to more significant uncer- 
tainties on the current constraints. We thus arrive at a 95% 
CL upper limit of ~ 2 — 3 nG on Bo. 



3.2 Constraints from the ionization degree 

Fig. [2] shows the constraints due to the observed ionization 
de gree, x e g, and its uncertainties. Following the analysis 
of iBolton et al.l (|201ll ). the IGM is assumed mostly ionized 
at z ~ 7. Note that for the present purposes the accurate 
value of the ionization degree is not crucial and even a 90% 
value would lead to very strong constraints independent 
of the specifically adopted CL. In our reference case (top 
panel), we find at z = 7 an ionization degree close to 1 for 
co-moving field strengths up to ~ 1.5 nG. For increasing field 
strengths, x e g drops rapidly, as the magnetic Jeans mass in- 
creases with Bq, significantly suppressing the production of 
ionizing photons. Thus considering all the uncertainties due 
to the cosmological parameters leads to Bo < 2nG. 

As already discussed for the optical depth, the uncer- 
tainties due to the reionization parameters are considerably 
higher (mid panel). However, cases with high clumping fac- 
tors, low escape fractions or l ow IGM temperatures and 
even very high esca pe fraction (|Srbinovskv fc Wvithel 120101 ; 
IMiniati et afll2004l ) are inconsistent even in the limit of weak 
magnetic fields. For the remaining cases, the ionization de- 
gree drops rapidly for Bo above a nG, and the most conser- 
vative upper limit corresponds again to Bo^ 2nG. 

As for the uncertainties in the parameters of the 
Schechter function bottom panel of Fig. [2] shows that some 
extreme values can be ruled out because inconsistent with 
the ionization degree at z = 7, even with negligible Bo. How- 
ever, the most important source of uncertainty for the upper 
limit on Bo are those related to My V . This forces our con- 
straint to B0S1 3nG. This constraint will be improved once 
the Schechter function is observed more accurately. Again, 
we checked that the uncertainties in the redshift evolution 
have no impact on our results. From the observed ionization 
degree, we thus find an overall constraint Bo ^2 — 3 nG. 

In Fig. [3]we show the ionization degree at different red- 
shifts as a function of the co-moving field strength. The plot 
illustrates the sensitivity of our constraint to the redshift 
at which the ionization degree is measured. In particular 
our constraint would be slightly weaker if based on data at 
z — 6, while it would considerably improve if a high ioniza- 
tion degree at z = 8 could be observationally established. 
At even higher redshifts, our model predicts the presence 
of a substantial neutral fraction in the IGM which could be 
translated into a further independent constraint. 



4 DISCUSSION AND OUTLOOK 

Using the combined constraints from the observed Thom- 
son optical depth, the ionization degree of the IGM at z = 7 
and the observed UV luminosity function at high redshifts, 
we have derived robust upper limits of 2 — 3 nG, virtually 
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independent of the CL, on the strength of primordial mag- 
netic fields. Previous work (e.g. Sethi & Subramanianll2005l ; 



iTashiro fc Sueivam a 2006; Schle icher et al.ll2008l . see also In- 
troduction) found similar limits but generally was based on 
the Thomson optical depth only, used less refined physical 
models, had a lower CL. However, the similarity of con- 
clusion also suggests that we are approaching the intrinsic 
limitation of the methods. 

Our results are mostly based on the lower cutoff of the 
luminosity function of galaxies contributing ionizing pho- 
tons, set by the magnetic Jeans mass (Eq. [7|). On the other 
hand, heating effects due to ambipolar diffusion are impor- 
tant for Bo close to 10 nG, which is ruled out by CMB ob- 
servations. Our constraints on the primordial magnetic field 
would be mostly strengthened by measurements of a high 
ionization degree at redshift zij8. This is simply because 
the earlier the IGM is reionized, the more important be- 
comes the contribution of ionizing photons from halos whose 
collapse is hindered by primordial magnetic fields. However, 
improvements beyond an order of magnitude are likely out of 
reach, as the formation of structure at the magnetic Jeans 
mass corresponding to 10~ 10 G (see Eq. [TJl is severely af- 
fected by other feed back processes (see, e.g., discussion in 
iBouwens et aljfeoill ). Other limiting factors in the deriva- 
tion of our constraints are mostly due to uncertainties in the 
UV luminosity function at high redshift and lack of infor- 
mation about the IGM ionization fraction at redshift higher 
than 7. There is reason to expect good progress in these ar- 
eas in the near future. For example, JWST[J will provide a 
unique opportunity to probe the IGM at and beyond red- 
shift 10, Lyman a emitters may provide additional informa- 
tion on the stren gth of the photoionizing background (e.g. 
lLatif et al.ll201ll ) , while the Planck satellit^f] will provide an 
improved measurement of the Thomson optical depth. These 
additional data may help to improve our understanding both 
with respect to reionization and the origin of magnetic fields. 
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